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For the formation of a black hole in the gravitational collapse of a massless scalar field, we calculate 
a critical exponent that governs the black hole angular momentum for slightly non-spherical initial 
data near the black hole threshold. We calculate the scaling law by second-order perturbation 
theory. We then use the numerical results of a previous first-order perturbative analysis to obtain 
the numerical value /j, ~ 0.76 for the angular momentum critical exponent. A quasi-periodic fine 
structure is superimposed on the overall power law. 



I. INTRODUCTION 



In the space of initial data for a self-gravitating system, the black hole threshold is the limit between initial data 
that eventually form a black hole and initial data that do not. More precisely, if one considers smooth one-parameter 
families of initial data such that for large values of the parameter p a black hole is formed, but not for small values, 
then there is a critical value p* (that can be found by bisection) where a black hole is first formed. Such data are 
called critical data. All numerical experiments are compatible with the picture that the critical data form a smooth 
y—i • hypersurface of codimension one, called the critical surface, in the infinite-dimensional space of initial data. 
t— I ' The evolution of any data sufficiently close to the black hole threshold shows what is now commonly known as 
critical phenomena in gravitational collapse. The solution evolved from near-critical data, with p ~ p*, approaches a 
universal intermediate attractor in a region of spacetime. This universal "critical solution" has a discrete or continuous 
self-similarity. Numerical experiments as well as perturbative calculations are compatible with the picture that it is 
C7^ an attractor inside the critical surface, with precisely one unstable perturbation mode pointing out of the critical 
^ | surface. 

If p > , and a black hole is formed, its mass scales with p as Q 

m*c(p- p .)i, (i) 

where 7 is the same for all families of initial data. However, scaling is not only a supercritical property. Any 
dimensionful variable of the system scales with p — p*. For example the maximum curvature of spacetime scales as 
(p* — p)~ 21 for p < (|]. The critical exponent 7 is directly related to the growth rate of the one unstable mode. 

The critical exponent 7 is independent of the initial data. It depends on the type of matter (or absence of matter, 
in the collapse of pure gravitational waves), but is independent of any dimensionful constants in the matter equations 
of motion. Critical phenomena were first observed in the spherically symmetric collapse of a scalar field ffl and the 
axisymmetric collapse of gravitational waves |s| . The structure of the intermediate attractor that plays the role of the 
critical point was clarified in pj and [H] . The critical exponent for the black hole mass was calculated by perturbation 
theory around the critical point in gTand more generally in Q. An independent critical exponent for the black hole 
electric charge was predicted in B and subsequently measured in The prediction of a critical exponent for the 
black hole's angular momentum required perturbation theory beyond spherical symmetry, and was first carried out 
in (n]Jll[] for perfect fluid collapse. A detailed review of all this can be found in p2f . 

Now two of us have calculated the non-spherical linear perturbations around the critical point in scalar field 
collapse [Q. In the present paper, we build on that work to quantitatively predict the scaling of the black hole 
angular momentum in critical scalar field collapse. The perfect fluid and scalar field cases have in common that the 
critical spacetime is spherically symmetric. Black hole angular momentum, which breaks spherical symmetry, can 
therefore be treated perturbatively. There are two important differences, however. 

The first difference is this. The critical point spacetime in perfect fluid collapse has, apart from spherical symmetry, 
another continuous symmetry, namely continuous self-similarity (CSS), also called a homothety. Roughly speaking, 
this means that it does not depend on t and r separately, but only through the dimcnsionless combination x = —r/t 
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(for a suitable choice of coordinates that we do not need to discuss here). The perturbations of this solution must 
then have the form Re [(— 1)~ f(x)] , where both the eigenvalue A and the mode function f{x) can be complex. In 
the perfect fluid case the A of the perturbation that dominates angular momentum - the I = 1 axial mode with the 
largest real part - is real. (The statement in |llf| that A is complex was wrong.) The black hole angular momentum 
in critical perfect fluid collapse scales as 

L ~ L (p - p*Y, (perfect fluid, CSS) (2) 

In contrast, the scalar field critical spacetime depends not only on x, but in a restricted way also on t: it is periodic 
in r = — ln(— t). (Again this holds only in suitable coordinates). This symmetry is called discrete self-similarity 
(DSS) [Q. Furthermore the critical exponent itself is complex. The imaginary part of the critical exponent and 
the independent period A in r of the background critical point solution then combine to give an overall power law 
modified by quasiperiodic behavior, as we shall see below. 

The second new difference between the fluid and the scalar field is that angular momentum can be treated in 
first-order perturbation theory around spherical symmetry for perfect fluid collapse, but not for scalar field collapse. 
We have to go to second-order perturbation theory to obtain a non-zero effect. Fortunately the relevant second-order 
perturbation degree of freedom does not obey a wave-like equation which would require independent free initial data, 
but is totally determined by its source, which is quadratic in first-order perturbations. This allows us to obtain the 
desired result with a calculation based on dimensional analysis and selection rules for the angular dependency of the 
perturbations. At the end of the calculation we insert the numerical value for a particular complex eigenvalue A into 
our result to obtain a numerical value for the angular momentum critical exponent. 

The paper is structured as follows. First we define perturbation theory around spherical symmetry to all orders and 
show, both formally and intuitively, why we have to go beyond first-order perturbations. Then we consider the second- 
order perturbations that give rise to angular momentum. In order to make our presentation more self-contained, we 
briefly review how the mass scaling law is derived. Then we derive the angular momentum scaling law for a spherical, 
DSS, background, and insert numbers for the scalar field case. 



II. FIRST-ORDER PERTURBATIONS AROUND SPHERICAL SYMMETRY 

In order to define higher order perturbation theory, we formally expand the metric and the scalar field as 

9»u(e) = + e fl W + e 2 g$ + 0(e 3 ), (3) 

0(e) = (o) + e0 (1) + e 2 (2) + 0(e 3 ). (4) 

In the following we use the shorthand u to denote both and 4>. If we write the field equations, that is, the Einstein 
equations and the scalar wave equation, formally as 

£(u)=0, u= (g^,^) =u(e), (5) 

and take formal derivatives with respect to e, we obtain as the leading orders 

£ (w (0) ) = 0, (6) 
C(u^) = 0, (7) 
£( M ( 2 >)=S( 2 W (1) ), (8) 

where £ is a linear derivative operator, and S is a quadratic derivative operator. 

We perturb around a spherically symmetric solution . In this paper we consider a double perturbation expansion, 
both around spherical symmetry, and around criticality. The small parameter e in the following always refers to 
deviations from spherical symmetry. (The small deviation from criticality will be parameterized by p — p„.) 

We begin by establishing that no black hole angular momentum can arise in first-order perturbation theory, for 
a spherically symmetric background solution with only scalar field matter. We begin by noting that the Kerr 
metric with mass M and angular momentum L in Boyer-Lindquist coordinates can be written as a perturbation of 
the Schwarzschild metric (for simplicity we assume that the angular momentum is in the z-direction): 

ff C° W dx v dt 2 + (l - ™\ 1 dr 2 + r 2 (dO 2 + sin 2 dtp 2 ) , (9) 

g^dx^dx" = -4- sin 2 9 dtdtp = 4-shi0^\ —Y w (6,<p)dtdcp, (10) 
p r r o9 V 3 
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where Yio is the spherical harmonic with I — 1 and m = 0. Note that this is the only first-order metric perturbation 
when we expand in powers of L. The other metric coefficients are changed only from 0(L 2 ) on. 

We adapt the gauge-invariant perturbation framework of Gerlach and Sengupta ]l5j , which is reviewed in more 
detail in |13| . In this framework, all four-dimensional tensor perturbations are covariantly decomposed into series 
of tensorial spherical harmonics, which carry all the angular dependence, with coefficients which depend on r and 
t. These harmonics are tensors derived from the Yj m and their derivatives, living on the two-spheres of spherical 
symmetry (coordinates 9 and (p, tensor indices a). The coefficients of those expansions are tensors living in the 
two-dimensional reduced manifold (coordinates r and t, tensor indices A). This decomposition is performed because 
around spherical symmetry different I, m components of the expansion decouple, so we can study each case separately. 
Each l,m linear perturbation decouple further into two parts: axial [with parity (— 1)' +1 ] and polar [with parity 

(-!)']■ . _ 
In this notation, the linear metric perturbation (nOh is written as 

gtydaP dx v = 2k A 1] dx A S 10a dx a , (11) 

where S±o a is the axial harmonic vector field with angular dependence I — 1, m = 0, formed by the general rule 
Sima — £a b Yi m :b- (Here : a is the covariant derivative on the unit two-sphere, and e a t, is the corresponding unit 
antisymmetric tensor.) The angular momentum perturbation L z we are interested in can therefore be characterized 
as the axial I = 1, m = 0, gauge-invariant metric perturbation, (m = 1 and m = —1 parameterize L x and L y .) By 
comparison with ([H]) we read off that 

k W dx A = ./±E— dt (12) 
V 3 r 

for I — 1 and m = 0. k^} is only partially gauge-invariant for I = 1, but its curl is gauge-invariant for all I > 0, 
and contains all the gauge-invariant information. We calculate 

^) = e^(r-^) {B = -^ 6 ^, (13) 

for I — 1 and m — 0. (In the notation of |D| these objects are simply called Ua and II because there we only work 
with first-order perturbations.) \A is the covariant derivative on the reduced manifold, and e AB the corresponding 
unit antisymmetric tensor. 

Let us compare this particular solution with the general linearized equation (^) it must obey. For I > 2, II^ 1 ) obeys 
a wave equation, but for Z = 1 it obeys a first-order differential equation that can be integrated trivially to yield 



r 4 



n (1) - 16ttT (1) - c = 0, (14) 



where is a scalar constructed from the axial gauge invariant matter perturbations, and c is an integration constant. 
Note that C degenerates from a derivative operator to an algebraic operator for I — 1 axial perturbations. In spite 
of its trivial appearance, equation ( |l4| ) is the linearized Einstein equation relating gauge-invariant axial I — 1 metric 
perturbations to their matter sources. 

If the spacetime is vacuum, the spherical background must be Schwarzschild, and the matter perturbation 
vanishes. The only axial physical perturbation is then the one that takes Schwarzschild infinitesimally into Kerr. We 
see that for this perturbation 




TW=0, c=-6J^-L. (15) 



(There is one other physical perturbation of Schwarzschild, / = polar, which changes the mass of the black hole.) 

If we demand regularity at the center r — 0, as we do here, we must have c = 0, and this is also true in the presence 
of matter. The scalar field and its perturbation (to all orders) are polar, simply because they are spacetime scalars. 
Therefore vanishes for scalar field matter and we have IT 1 ) = 0. (To the next order, polar and axial perturbations 
do mix, so that to second order T and hence II do not vanish. This will be discussed in the next section.) 

While this is the complete argument for the absence of black hole angular momentum in the first-order perturbation 
calculation, we have also found two intuitive ones. Angular momentum is present in the first-order perturbations of 
spherical fluid collapse, but not of scalar field collapse, because the fluid is made up of individual particles with 
a rest mass which can go round while the entire configuration remains axisymmetric. The angular momentum is 
proportional to the tangential velocity component u v . The scalar field has no such particles, and we need to make 
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spoke-like structures in the field (thus breaking axisymmetry) and then make these go round. Intuitively, we need 
two powers of e to do this: one to make spokes, and another to make them go round. This argument is backed up by 
the observation that the conserved angular momentum (Noether charge) of a scalar field on a Minkowski background, 
for example the z-component, is simply 

L l, Minkowski = J M,<p d 3 x = 0(e 2 ). (16) 

One can easily see that this is again quadratic in deviations from spherical symmetry. 

A second intuitive argument for the absence of first-order angular momentum perturbations in the scalar field 
comes precisely from its partial equivalence with a perfect fluid (with p — p): the 4- velocity of that pseudo-fluid is 
the normalized gradient of the scalar field, and therefore irrotational. 



III. SECOND-ORDER PERTURBATIONS 



In second-order perturbation theory, we can define second-order versions of the gauge-invariant perturbations. 
(Equivalently, we could fix the gauge separately at each order.) Independent of the detailed field equations, the term 
in the n-th order equations that is linear in the n-th order perturbations is the same at all orders, namely C. For 
I = 1 axial perturbations, equation (||) therefore takes the simple algebraic form 

r 4 n( 2 > - 16nT^ - c = S« (u^u^) (17) 

where Il( 2 ) and T< 2 ) are just IT 1 ' and T^ 1 ' with replacing u^ 2 \ Just as for the first-order perturbations, 
vanishes identically for scalar field matter, and the integration constant c vanishes if we consider only solutions with 
a regular center r = 0. 

But now there is also the source S^ 2 \ which is quadratic in first order perturbations It does not vanish in 

general, and generates a non- vanishing II^ 2 ^ . Therefore it can introduce angular momentum. If and when a black hole 
is formed, it must settle down to the Kerr solution at late times. must then approach a constant at late times 
outside the horizon, thus transforming a Schwarzschild into a Kerr black hole. Inside the horizon, where the singularity 
forms, perturbation theory around a regular center will break down, but that does not affect our calculations. 

Finding the expression for would require more effort than writing down the left-hand side, but fortunately we 
are interested only in scaling arguments, not in the detailed behavior of ±l( 2 ) as a function of r and t. We note that II^ 2 ) 
belongs to the axial sector, and that we are only interested in its I = 1, m = component, as that is the one connected 
to black hole angular momentum. All other second-order perturbations (except a mass perturbation, as mentioned 
above) must eventually be radiated away as the black hole settles down to the Kerr solution. We can therefore restrict 
attention to terms which are quadratic in first-order perturbations, axial, and have angular dependence I = 1, m = 0. 
The restriction to m = is equivalent to the restriction to the z-componcnt, L Zl of angular momentum. We make it 
only for simplicity of presentation. L x and L y are related to complex linear combinations of m = 1 and m = — 1. 

Let us denote the two factors u' 1 ' by u' and u" . They must have m — ml + m" = 0. Wc do not require that the 
number I' + I" is odd, nor that of v! and u" one is polar and the other axial. An example for the mixing of polar and 
axial perturbations to quadratic order is 



Yim Yj_ m:a = -im \ —— S Wa + other terms. (18) 

V 1D7T 

for any I and m. Note that any m' + m" = gives rise to the axial vector field Sioa that characterizes angular 
momentum in the z-direction, except m! = m" = 0. That ml — m" = is excluded is also clear on physical grounds, 
as we have argued above that an axisymmetric scalar field configuration cannot have angular momentum. In the 
simplest case I = 1 the complete expression is 



Y n Y 1 _ 1:a = -iJ^S Wa + -^=Y 20 .. a . (19) 

V l07T 4V57T 

In words: two polar I — 1 perturbations combine to give an axial I = 1 perturbation, as well as the expected I = 2 
polar perturbation. 
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IV. REVIEW OF MASS SCALING IN CRITICAL COLLAPSE 



We have clarified the role of second-order perturbation theory for calculating angular momentum- like quantities in 
almost spherical scalar field collapse. In the following we specialize to a particular background solution, the so-called 
critical solution. The derivation of the angular momentum scaling law is a continuation of the previous derivations 
of scaling laws for the mass QJt]] , electric charge B and angular momentum . Therefore we do not try here to be 
fully self-contained, but rather remind the reader of the general ideas underlying these calculations. A more detailed 
review is contained in fll2| |. 

Critical collapse is dominated by a single solution which has the two crucial properties of being self-similar (CSS 
or DSS) and of having precisely one growing perturbation mode. This solution is best given in coordinates x and t, 
which are (with some simplification) x = —r/t and r = — ln(— t). The important properties of these coordinates are 
that both are dimensionless, that x is invariant under rescaling of space and time, and that r, being the logarithm of 
a scale, changes only by an additive constant. The critical scalar field solution, in coordinates t, x, 9, tp is then of the 
form 

g liv = lle~ 2r g lliv {x,f), <j> = ()>{x,t), (20) 

where the tilde over the argument r indicates that the conformal metric g^ and scalar field (f> depend on r only 
periodically, with a period A. Note that g^v, <f>, and the coordinates r, x, 9, tp are conveniently thought of as dimen- 
sionless, with the arbitrary constant scale Iq the only dimcnsionful quantity. Note also that by construction, Iq always 
comes together with e~ r . From this it follows that in the self-similar spacetime 

any masslike quantity ~ loe T - (21) 

Now consider the one growing linear perturbation mode. As the background is periodic in r, it must be of the form 

growing mode ~ (p - p*)e A ° T f (x, f). (22) 

Ao is real and positive, so that the perturbation grows as we approach smaller scales (t — > 0_, r — > 0, t — * oo). As 
by assumption there is only one such mode, and the background is real, Ao and / must be real, but in general Ai 
and fi exist in complex conjugate pairs. The overall amplitude of the growing mode depends on the initial data in a 
complicated way, but to leading order it must be proportional to p — p* , because for p = p* the critical solution lives 
forever, by definition of being the critical solution, and so the growing mode cannot be present. 
Now consider the spacelike hypersurface r = r*, with the value of t* defined by 

growing mode ~ (p — p»)e A ° r * ~ some fiducial amplitude. (23) 

For r* this gives 

A t* = - ln(p - p*) + const. (24) 

At a later stage, we can no longer approximate the spacetime as self-similar plus a perturbation, but the Cauchy data 
at r = r* are independent of the initial data - the decaying perturbations have all decayed, and the one growing 
perturbation has reached its fiducial amplitude - up to an overall scale, which must be given by @. Therefore we 
have for the black hole mass 

M~ he- T » - (p-p.)*, (25) 
so that we have found the law ([!]) with 7 = 1 / Ao . 

V. ANGULAR MOMENTUM SCALING IN DSS CRITICAL COLLAPSE 

At the fiducial time r = t* , we have for the angular momentum of what will become the black hole 

L z ~ (r*) ~ l 2 e- 2T *Reu'{T*) Reu"(r*) ~ l%e- 2r *Re[C /'(r*)e A ' T «] Re[C"7"(n)e A " T «] (26) 

Here A' and A" are the complex exponents of the most slowly decaying modes of two separate first-order perturbations 
which are compatible with the angular dependence selection rules, and /' and /" arc the corresponding complex mode 
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functions, which are periodic in r with period A. C and C" are complex constants that depend on the family of 
initial data. (They also depend on p, but here we only take their leading order, which is constant. We take the real 
part of v! and u" separately because the complex notation is only a shorthand for sines and cosines.) In order to 
simplify we have made the jump from the x-dependent quantity to the simple number L in the first equality, 
and we have therefore suppressed the x-dependence of /' and /". Note that the factor Iq e~ 2r * in the second equality 
appears by dimensional analysis - L has dimensions (length) 2 . 

In [T^ | we found that the most slowly decaying perturbation with / 7^ is the polar I = 2 perturbation, which has 
A ~ —0.06 x (1/A) + 0.30 x (2ni/A) ~ —0.017 + 0.55i. Therefore L z will be dominated by contributions where both 
v! and u" are I = 2 polar. This does not mean, however, that u' = u" . will contain dominant contributions 

from both ml — —ml' = 1 and m! = —m" = 2. Furthermore C and C" are independent complex constants for each 
of these two cases. (We have chosen our complex notation so that C can take arbitrary complex values, but A and 
its complex conjugate are not counted as independent.) The function L z (p — p*) is therefore parameterized by four 
independent complex constants, namely the constants C for I — 2, polar, with m = — 1, 1, —2, 2. 

Inserting the value (|||) for t* , and putting back the vector character of L, and the periodic nature of the background, 
we can summarize our result as 

L^LoKp-pOKp-P.)". (DSS) (27) 

where Lo[ln(p — p*)] is a quasiperiodic function that depends on the family of initial data, and fi is a universal critical 
exponent that is given by 

fi = Rc 2 ~\ - A = (2 - 2ReA')7- (28) 
Ao 

Note that A' has negative real part, so that we have p > 2j. With 7 ~ 0.374 and ReA' ~ —0.017 we predict a critical 
exponent /1 ~ 0.76, which is barely larger than 27. 

The Fourier spectrum of Lo[ m (p — P*)\ with respect to its formal argument contains the angular frequencies 

N^, iV^±2 7 ImA' (29) 

for integer N. We could be more precise by writing down the general form of Lo[ln(p — p*)] as a sum involving eight 
universal (but from the present calculation, unknown) periodic functions and 24 real constants depending on the 
family of initial data, but that would not be very helpful. Nevertheless, in numerical collapse simulations it should 
be possible to spot not only the overall power law, but also a fine structure with the fundamental angular frequency 
27T7/A ~ 0.683 and the offset 271mA' ~ 0.41 in the Fourier transform of Lq with respect to ln(p — p»). 
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